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A LIGHTING STRATEGY FOR LUNAR ORBITER MISSION DESIGN 


By Friedrich O. Huck 
Langley Research Center 

SUMMARY 

An analysis is presented which relates quality of lunar photographs to Lunar 
Orbiter trajectory. Furthermore, since several parameters which affect picture quality 
are uncertain, a strategy is developed which optimizes confidence levels in detecting 
specific lunar surface detail. 

This strategy is illustrated by examining Lunar Orbiter's capability to reveal small 
obstacles in the form of right- circular cones which have lunar reflective properties. 
Results show that the smallest lunar detail can be detected at high phase angles (angle 
between incident light and camera) and long exposure times; however, the effects of 
variances in lunar reflectance, vehicle trajectory, and image smear result in higher 
levels of confidence in the detection of larger detail at lower phase angles and shorter 
film exposures. 


INTRODUCTION 

The prime objective of the Lunar Orbiter project is to obtain topographic informa- 
tion regarding various lunar areas to assess their suitability for use as Apollo landing 
sites. Since photographic quality in revealing small obstacles which would be hazardous 
to the Lunar Module varies strongly with vehicle trajectory and film exposure, a criterion 
is needed for the selection of these parameters to optimize detail detectability. 

Several measures of photographic quality have been established to define Lunar 
Orbiter's capability to aid in the Apollo program. This paper presents an analysis 
which relates these measures of picture quality to viewing geometry and camera shutter 
speeds. Furthermore, this paper considers the uncertainties which exist in lunar reflec- 
tive properties, vehicle trajectory, and camera system operation by developing a 
lighting strategy for securing maximum levels of confidence in the detection of specific 
lunar features as a function of viewing geometry and camera shutter speeds. 

An image signal-to-noise ratio of lunar features is selected as an appropriate 
assessment of picture quality because it includes the effects of lunar reflection, viewing 
geometry, and camera characteristics and because it is a measure of the information 


content in the photographic image. Because of the high cost of obtaining closeup photo- 
graphs of the lunar surface, it is desirable that these photographs contain maximum 
information about small detail. Thus, a convenient measure of Lunar Orbiter's capa- 
bility to reveal small depressions and proturberances has been established in terms of 
the detectability of right- circular cones which are near the limits of the camera system 
resolution. This measure of picture quality is introduced in this paper to illustrate the 
proposed lighting strategy. 

Because the detection of small detail is of primary interest, results of this analysis 
are given for the high- resolution camera system and an expected vehicle perilune of 
46 kilometers. Under these conditions, Lunar Orbiter's vertical viewing geometry can 
be expressed as a function of phase angle alone, where phase angle is defined as the angle 
between incident light and camera axis. And, because lunar targets near the limits of 
the system resolution are considered, the system may be linearized without introducing 
significant errors. 

While these factors simplify computations, the basic analysis is not restricted to 
them, but represents a general approach to the prediction of picture quality and optimum 
viewing geometry for the photography of moon and planets by a spacecraft camera. 

SYMBOLS 

A area: on lunar surface, meter 2; on film image, millimeter ^ 

B brightness, meter-candles 

d cone diameter: on lunar surface, meters; on film image, millimeters 

E exposure, meter-candle- seconds 

f n focal ratio or "f-number" 

F spread function 

g phase angle, degrees 

G modulation transfer function (MTF) 

h height, kilometers 

k fundamental spatial frequency component of cone model, lines/millimeter 
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K 


solar constant, meter- candles 


L linear operator 

m,n dummy variables 

N root-mean- square (rms) noise 

O lunar object 

p probability density function 

s smear, la value of Gaussian smear distribution: on lunar surface, meters; 

on film image, millimeters 

5 signal 

S/N signal-to-rms-noise ratio (hereinafter referred to as signal-to-noise ratio) 

t exposure time, sec 

T film transmission 

u,v spatial frequency, lines per millimeter 

x,y,z spatial dimensions: on lunar surface, meters; on film image, millimeters 

/3 lunar slope 

y confidence level 

6 difference between actual mean value of minimum detectable cone diameter 

and its estimate 

r nonlinear operator 

9 angle, degrees 

A vector parameter set representing spatial dimension or spatial frequency 
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H actual mean value 

If vector parameter set containing all variables which affect lighting strategy 

p lunar albedo 

a standard or rms deviation of a Gaussian distribution 

t static lens transmission 

<f>( ) photometric function 

< > expected value 

Subscripts: 

a aperture 

b background 

c film readout and communication subsystem 

D film density 

d cone diameter 

e effective 

f film 

i image 

l lens 

max maximum 

min minimum 

p peak 


4 



q dark 

r bright 

s smear 

T film transmission 


ANALYSIS 

Object Detectability Threshold 

Signal analysis .- The most general description of the purpose of a camera is per- 
haps that it is to give information about an object. If it is required to extract all the 
information that is implicit in the resulting image regardless of the possible complexity 
of the interpretation, the quality of an image is most appropriately assessed by its 
signal-to- noise ratio (ref. 1). Accordingly, an image signal-to-noise ratio of lunar 
features is derived to predict picture quality. Since the detection of small objects is of 
primary interest, it is convenient to give this ratio for signals near the limits of system 
resolution. For this purpose, a simplified block diagram of the Lunar Orbiter photo- 
graphic system is given in figure 1. The properties of lunar reflection (refs. 2 to 5) and 
the assessment of optical images in the presence of signal degradation and noise (refs. 1, 
6, 7, and 8) are discussed in this analysis only to the extent to which the proposed 
lighting strategy and hence photographic mission design is affected. 

All system elements, except film response to exposure, can be considered linear 
(ref. 9) ; film conversion of exposure to transmission can be treated approximately as 
linear only if very low contrast objects are considered. It is convenient here to define 
A as a vector parameter set representing spatial dimension or frequency variables and 
t as a vector parameter set containing all variables which affect lighting strategy, and, 
furthermore, to define L as a linear operator and r as a nonlinear operator. 

For an object description 0(A ;/3) with spatial dimensions or frequencies A and 
slope /3, the photometric function of the lunar surface yields a brightness variation 
which can be expressed as (ref. 2) 

AB(5C;j3,g) =Kp A<fr(Y;0,g) (1) 

where 

K solar constant at moon surface 

p full-moon albedo or reflectance coefficient for normal incidence and emission 
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*0,g) 


surface photometric function normalized such that $(0,0) = 1 


This brightness variation is a function of viewing geometry. For the Lunar Orbiter's 
vertical photography, viewing geometry may be expressed as a function of phase angle g 
as defined within a plane in figure 2(a). Brightness contrast as a function of a more 
general viewing geometry is given in reference 10. 

Exposure variations at the spacecraft camera focal plane are related to brightness 
by (ref. 2) 


AE(r i; /3,g,t) = AB(r;0,g) (2) 

4f n 

where 

t exposure time 

f n camera lens f-number 

r static optical transmittance 

Lunar object dimensions X and film image dimensions xTj are related by camera 


lens focal length and vehicle altitude over the object. Exposure variations may, there- 
fore, for a given lens be written as a function of vehicle height h as AE(r;0,g,t,h). 

Because of degrading effects due to image smear s and lens-film characteristics, 
the actual film exposure variation is 

AEf(r ; ?) = L z>f {r)L s (x;s){AE(r ; /3,g,t,h)j 

( 3 ) 

where if = jS,g,t,h,s. This equation may be rewritten as 


AE £ (r ; ?) = q jfjS (^;s){AE(r;ftg,t,hj} 

( 4 ) 

where 


L Z,f,s( X ’ s ) = L s(^; s ) 

( 5 ) 


If the vector parameter set X is defined as a spatial dimension set (x,y), then 
the operation of L, f (x,y;s) on the input signal AE(x,y:/3,g,t,h) becomes the con- 

6,1, S 

volution integral 

AEf.(x,y;lf) = J J ^ F^ fjS (x-x , ,y-y';s)AE(x’,y , ;/3,g,t,h)dx' dy' (6) 
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where £ g (x,y;s) is the combined lens-film and image smear spread function. If, 
however, the vector parameter set X is defined as a spatial frequency set (u,v), then 
£ g (u,v;s) indicates multiplication, and equation (4) can be written as 


AE f (u,v;|) = G l>f)S ( u ,v;s) AE(u,v;/3,g,t,h) 


(7) 


where G. f (u,v;s) is the combined lens-film and image smear modulation transfer 
function (MTF). Spread function and MTF are related by the Fourier transform pair 


G(u,v) - f " f " F(x,y) e " i2ff ( ux+v y) dy (8) 

oo 


F(x,y) = f ~ f " G(u,v) e 12 ’ r <»*«'y>du dv 

^ -oo 


(9) 


Once film exposure is determined, the resulting film transmission can be evaluated 
from the film characteristic curve; this is depicted by using the nonlinear operator Tf, 
so that 

AT(jT;F)= r f (?){AE(r ; ?)} (10) 

The signal passes then from the film to the film readout and communications link 
which are combined here into one block. The output signal is related to film transmis- 
sion variations by the linear operator L c (r), as given by 

s(r ; |')=L c (r){AT(r ; ?)} (n) 

which, again, indicates convolution in the spatial dimension domain and multiplication in 

the spatial frequency domain. 

Combining equations (4), (10), and (11) and being careful not to disturb the correct 
order of operation with respect to the nonlinear operator gives the output signal in terms 
of the input exposure as follows: 

S()T;?)= L c (r)rf(?)L i>f>s (r;s){AE(r;ftg,t,h)} (12) 

When scenes near the system resolution limit, yielding images of low contrast, or scenes 
which have low contrast are considered, the image distribution can be expressed as 


s(x,y;? ) = J J AE(u,v;0,g,t,h) c («» v J s )^^)e l2ir ^ ux+vy W dv (13) 
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where the linear operators L C (X) and f s (^;s) are replaced by the total system 
MTF Gj j s c (u,v;s), and the nonlinear operator is replaced by the linearized 

film response ^^(^b)* K ma y n °t e( i that film response now becomes a function of 
the background exposure E b . 

Noise analysis .- All system noise sources are assumed to be "white," statistically 
independent, and represented by a Gaussian process. The total rms system noise is 
grouped for the present purpose into exposure- dependent film granularity noise and 
exposure-independent communications noise. A detailed discussion of all system noise 
sources is given in reference 9. 

Film granularity noise is measured as the standard or rms deviation about the 
mean density and is commonly designated by ctj-j. This rms value depends on the actual 
film graininess and the film scanning aperture area (ref. 11, p. 53) which is treated 
herein, in harmony with the discussion in the previous section, as part of the communica- 
tion link. The rms fluctuations in film transmittance Vrj, may be obtained from 
by the relation (ref. 12) 

a T ( E b) = 2.3T(E b )a D (E b ) (14) 

where T is the average film transmission. 

Noise in the communication link can be expressed in terms of an equivalent value 
a c at the system output, and it can be root sum squared with film transmission noise to 
yield a total rms noise as given by 

a ( E b)= [V( E b) + CT c 2 ] V2 (15) 

But noise content in an image depends not only on system noise sources and MTF, 
it depends also on the signal image area. This is simply because the larger the image 
area is, the more samples are obtained by the readout aperture as it scans across the 
film. The number of independent samples n that may be obtained from an image area 
Ai by an aperture area A a is A^/Aa- According to the Central Limit Theorem 
(ref. 13), the mean of a sample of size n has a standard deviation of cr/Vn if cr is 
the population deviation. The standard deviation from the mean value of an image area 
larger than the scanning aperture area is, therefore, 



Signal-to-nois e rat io .- Combining equations (13) and (16) gives an image signal-to- 
noise ratio for targets near the system detectability threshold as follows: 
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§(*,7i?) = ^£(Eb )|g J_” AE(u,v; ft g,t,h) G M , SjC (u,v; S )e i2 ”< ra+ ^>du dv (17) 


with exposure and frequency dependent system characteristics separated. 

Object detection thresholds depend, of course, not only on image signal-to-noise 
ratios but also on the process of image interpretation, and become, therefore, a quantity 
of some arbitrariness. Detection thresholds have, nevertheless, been established. 

From an experimental study on the detection and identification of geometrical shapes, 
such as right- circular cones, in a simulated lunar environment (ref. 5) the following 
rule evolved: An object may be detected if its film image peak-to-peak signal-to-noise 
ratio is above 3 and identified if its size is about three times as large as needed for its 
detection. This signal-to-noise ratio detection threshold may be slightly optimistic in 
applications to actual lunar photographs, because the photo interpreters had been told 
the exact location of the test targets in the photographs. In evaluating photographic 
image quality, Schade (ref. 14) found that small image elements can be detected on a 
uniform background if their signal-to-noise ratio is about 3 to 4. Another study on 
imaging systems and human vision by Rose (ref. 15) indicates a signal-to-noise ratio 
of 5 as the lower limit for object detectability through visual observation. 

However, the importance of calculating object detectability does not lie solely in 
the absolute magnitudes of the results. It is the variation of these results as a function 
of mission controllable variables which allows optimization of detail detectability by 
mission design. Because the Lunar Orbiter Project Statement of Work specifies as a 
mission goal that a cone of 2-meter diameter and a base-to-height ratio of 4:1 shall 
result in an image peak-to-peak signal-to-noise ratio of 3, this ratio is used as effective 
cone detectability threshold. 


Lighting Strategy 

Whereas equation (17) presents a fundamental relationship for calculating target 
signal-to-noise ratio, uncertainty of its results prevails due to variations in film expo- 
sure, phase angle, vehicle height, image smear, and other factors. Film exposure varia- 
tions are expected because of uncertainties in lunar reflective characteristics and photo- 
graphic subsystem properties; phase angle and vehicle height variations are expected 
due to trajectory uncertainties, and image smear variations are expected due to uncer- 
tainties in velocity- over- height sensor errors and other image motions (ref. 16). Of 
these parameters, optimum lighting is most critically affected by variations in exposure 
and phase angle. Image smear influences lighting primarily through the selection of 
exposure speed; camera height has no effect, but changes in height during exposure will 
cause an image smear known as the zoom effect. 
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A well-founded basis for selecting a lighting strategy must weigh these uncer- 
tainties. The strategy proposed here optimizes confidence levels in the detection of 
specific surface features. To accomplish this, an estimate of the actual, but a priori 
unknown, mean value of the minimum detectable target feature is made and used to place 
confidence limits on the detection of specified target features as a function of phase angle 
and camera shutter speed. The proposed lighting strategy is, for convenience, dis- 
cussed in terms of minimum detectable cone diameters. But the strategy is applicable 
equally well to any other type of target. 

Expected va lue of minimu m det ec table con e diameter.- As stated in the introduc- 
tion, the primary objective of the Lunar Orbiter project is to gather information on the 
detailed surface structure of various lunar areas. In particular, several (on the average 
about 10) potential Apollo landing sites are to be photographed during a mission. And 
several photographs (on the average about 16 frames, each containing a high- and 
medium- resolution picture) are taken of each site. 


In designing a phase angle (trajectory) and camera shutter speed for taking n 
photographs of a site, it is desirable that the average value of the detectable cone diam- 
eters for these n pictures be minimized. This average value is defined as 
n 


<^d> = H 2 M di’ where 
i=l 


is the a priori unknown minimum detectable cone diameter 


for the ith picture. The best available estimator of this quantity is the expected value of 
d m in(?)» which may be written as (refs. 13 and 17) 


^min/* ~ H ^ ^min(^) p( 


(18) 


This is a best estimate of <Md> * n the sense that <(d m i n ^> can be expected to 
occur most frequently. In general, (d m j n ) * < lmin( < (^))> where indicates the 

most likely values of J; these quantities are equal only in a linear system. In the case 
at hand, the lunar photometric value varies nonlinearly with phase angle, and the film 
transmission varies nonlinearly with exposure. 

From the viewpoint of designing spacecraft trajectory and camera shutter speed, 
it is not of primary interest to predict the actual mean value of the minimum detectable 
cone diameter but to predict an optimum trajectory and shutter speed. For this purpose, 
equation (18) may, for numerical computations, be stated in a simplified form as a func- 
tion of only those parameters which affect lighting geometry most critically, namely, 
film exposure and phase angle. The influence of image smear on viewing geometry and 
shutter speed selection can be observed by repeating these calculations for probable 
amounts of smear. Even though it is still convenient to talk about an expected value of 
the minimum detectable cone diameter, it should be kept in mind, however, that the 
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calculated mean values and confidence limits are actually still a function of image smear, 
vehicle height, and other parameters which affect object detection. 

With these considerations, equation (18) becomes 

( d min> = f f d min( E >g) P(E,g)dE dg (19) 

' *^-. 00^-00 

Confidence levels in detection of a specified cone diameter .- A fundamental goal of 
this paper may be stated as follows: If ^p^) is the actual, a priori unknown, mean 
value of the minimum detectable cone diameter of a finite number of pictures taken of a 
lunar area, and (d m i n ) is its estimated value, what confidence level (Conf.) y can be 
placed on the difference (p^) - ( d min) being equal to or less than some specified 
value 6? This statement can be expressed in the form (refs. 13 and 17) 

Conf.^(p<j) " ( d min) = s) = Y (20) 

If the value 6 is chosen in such a manner that (d m i n ) + 6 = d, where d is a 
specified cone diameter, equation (20) may be rewritten as 

Conf .(<Pd) = d ) = r ( 21 ) 

and read as the confidence level y that the mean value of the actual minimum detectable 
cone diameter is less than or equal to some specified diameter. And the confidence level 
of this event is, in turn, equal to the probability (Prob.) that the estimated values of all 
parameters which affect (V-d) are favorable to this event - that is, 

Conf. ((Md> s d) = Prob. (<? mln > s T s <? max >) = < 22 > 

where the range of t° (^max) includes all the estimated values of £ for 

which (V-d) = d * 

Again if only film exposure and phase angle variations are considered, equation (22) 
is rewritten as 


Conf./^p^ = d^ - Prob.^(E m i n ) = E § (E m ax); (gmin) = g = (gmax)j 

C (Smax) 

(gmin) 


C ^rnax) , 

\ . p(E,g)dE dg 

(Emin) 


(23) 
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If 


In essence, this equation is a summation of all the probabilities of those combinations of 
exposure and phase angle which yield an average detectable cone diameter less than a 
specified value d. 

AN APPLICATION OF THE LIGHTING STRATEGY 

The effect of the proposed lighting strategy on the selection of an optimum viewing 
geometry for detecting small cones is illustrated in this section. For this purpose, an 
approximate expression is derived which enables rapid estimates of the minimum 
detectable cone diameter as a function of mission variables. It is convenient to perform 
all calculations at the airborne film output. Signal degradation resulting from the com- 
munication link may be neglected without introducing significant errors in the general 
conclusions drawn from these calculations. However, before this formula can be 
developed, some properties of the lunar surface, the target signature, and the Lunar 
Orbiter photographic subsystem must be introduced. Computer programs have been 
written to calculate the detectability of cones and other targets with more detail and 
rigor (refs. 9 and 18); however, their results are in close agreement with those obtained 
here. And, while it is believed meaningful to present major system characteristics to 
illustrate computations, it does not aid the purpose of this paper to discuss all the system 
details that have been considered in the computer programs. 

Lunar, Cone, and System Properties 

Viewing geome try an d lunar r efle ctivity .- For the vertical photography of the 
Lunar Orbiter, viewing geometry and lunar slope can be defined in a plane as shown in 
figure 2(a). Lunar reflective characteristics and film exposure are related for this 
geometry by (ref. 2) 

E(ftg,t) = 4>(ftg) (24) 

«n 2 

The nominal values of the given parameters, as well as their tolerances, are given in 
table I. Curves of the lunar photometric function based on work by Hapke (ref. 3) and 
Willingham (ref. 4) are shown in figure 3. 

Exposure and phase angle probability variations are assumed to follow a Gaussian 
distribution. This is suggested by the fact that these variations depend on a number of 
independent occurrences of which many have already been assumed to follow this dis- 
tribution (ref. 16). 

Uncertainties in the parameters which determine film exposure are expressed as 
exposure variations and are root sum squared to yield an exposure tolerance. The 
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probability density function p(E,g), as given in equations (19) and (23), can be repre- 
sented, in general, by the product of p(E/g) and p(g). Though exposure is a function 
of phase angle, exposure uncertainty is not; hence, the form of the conditional proba- 
bility density function p(E/g) does not change with phase angle. Since exposure devia- 
tion is given as a percentage of its mean value, its probability density function is unsym- 
metric and changes with mean value. It is, therefore, convenient to express this density 
function in terms of log exposure in order that its shape be symmetric and not vary with 
changes in mean log exposure. A discrete approximation of the log exposure and phase 
angle probability density function is given in figure 4. 

Lunar c one model .- The Lunar Orbiter Project Statement of Work has described 
among other goals a 26.6°, 2-meter-diameter cone as a measure of small obstacle detect 
ability. An approximate model of this target is used in this paper to calculate its mini- 
mum detectable diameter. This model is depicted in figure 2. 

The value of the photometric function $ and, hence, lunar reflection depends on 
phase angle g and lunar slope /3. These angles are defined in figure 2(a) in a plane 
parallel to the incident sunlight and the light reflected toward the camera. Away from 
this plane, the effective lunar slope - that is, as lunar slope affects the photometric 
function - becomes a function of the angle 6 which is defined in figure 2(b). An expres 
sion which relates the effective lunar slope to the angle 9 is derived in the appendix. 

Since cone sizes near the limits of system resolution are of interest, the value of 
target brightness contrast with background can be taken as an average value (ref. 10). 

The choice of average values to represent a more complex cone surface brightness dis- 
tribution is not unrealistic because small details will be largely obscured. The average 
photometric value is calculated here separately over those regions of the image which 
are brighter than the background # r and those regions which are darker than the back- 
ground l>q. Results are given in figure 5. The estimate of $q becomes increas- 
ingly conservative as phase angles rise above 67° because effects of shadow are not fully 
considered. 

To obtain maximum signal-to-noise ratio for a deterministic signal in the presence 
of white noise, a filter should be matched to the signal (ref. 19). Correspondingly, an 
optimum scanning aperture would be matched to the cone model's film image. This is 
easier in theory than in practice. Hence, a mathematical aperture area is used which 
is representative of an actual densitometer; namely, two neighboring circular areas 
whose diameters are equal to the radius of the actual cone and whose centers are located 
along sunline. This is indicated in figure 2(c). It is convenient to consider this con- 
figuration as a cone model and assign the average photometric value for the bright and 
dark side of the cone to these areas. 
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The scaling factor between lunar object dimension and airborne film image for an 
expected vehicle perilune of 46 km and a lens focal length of 610 mm is 0.0132 mm on 
the film per meter on the lunar surface. 

Photographi c subsystem characteristics and image smear.- It was shown in 
deriving equation (17) that for a small signal analysis the photographic subsystem of the 
Lunar Orbiter can be grouped into a spatial frequency and an exposure-dependent trans- 
fer function and that the transfer function for smear can be combined with the 
frequency-dependent system characteristics. 

Frequency- dependent system characteristics: Signal degradation due to 
frequency-dependent system characteristics depends primarily on the lens and film 
modulation transfer function which is shown in figure 6(a). Noise generated in the com- 
munication link is nevertheless included. 

Image smear is assumed random since it is caused by a number of independent 
motions (ref. 16). Its effect can be predicted by use of the modulation transfer function 
(ref. 8) as given in figure 6(b). The curves shown are labeled by their standard devia- 
tion referred to the airborne film and the lunar surface for a vehicle perilune of 46 km. 
The combined lens and smear modulation transfer functions are shown in figure 6(c). 

Exposure-dependent system characteristics: Film transmission and granularity 
noise are exposure dependent. Nominal characteristics of film transmission plotted 
against exposure are shown in figure 7(a). Film transmission noise a-p and a root- 
sum-square (rss) value of this transmission noise and communication noise cr c (corre- 
sponding to an expected peak-to-peak signal-to-noise ratio of 30 dB) is plotted in 
figure 7(b). 

For small signal calculations, a film-sensitivity-to-rms-noise ratio is given in 
figure 7(c). Film sensitivity AT/AE was calculated for an exposure change of 
0.1 meter-candle-second because images near the limits of system resolution exhibit 
approximately this film exposure variation. 

Scanning aperture: Scanning aperture size, which determines the amount of film 

transmission noise due to granularity, also determines the number of samples obtained 
from a given image area. While the Lunar Orbiter photographic subsystem uses a 
Gaussian spot, Eastman Kodak Co. (ref. 16) measures granularity using a circular 
aperture of 6.48-micron diameter, and it is also convenient to use this circular aperture 
here. It has been shown in appendix D of reference 9 that granularity noise measure- 
ments yield equal results for these two apertures if the diameter of the circular aperture 
is 4a and a is the standard deviation of the Gaussian spot. 
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Calculations 


Minimum detectable cone diameter .- An approximate expression is now derived 
from equation (17) to enable rapid estimates of minimum detectable cone diameters as a 
function of mission variables. It is convenient, at first, to consider only the signal term. 

The signal component of equation (17) is 

S(x,y;g,t,s) = (E b ) y ^ J ^ AE(u,v;g,t) dv (25) 


where the film readout and communcations system MTF are omitted. This expression 
is no longer a function of the variables 0 and h since these variables have been 
specified previously by selecting a 26.6° cone target and an expected vehicle perilune of 
46 km. The brightness distribution of the target model was represented by two neigh- 
boring areas with an average photometric value of the bright and dark side of the cone. 
Because equation (25) is linear, the rule of superposition allows the two areas to be 
handled as separate inputs to the camera system and to be combined again at the output. 
Thus, the peak output signal of each area, S r p(g,t,s) and Sqp(g,t,s), can be determined 
and later combined to yield the desired peak-to-peak output signal 


Sp_ p (g,t,s) = S rp (g,t,s) + Sqp(g,t,s) 


It is mathematically convenient to locate the X,Y coordinate axis in such a man- 
ner that the signal area is symmetric around it. The signal content of, say, AE r (x,g;g,t) 
is then averaged over an area bounded by -x 0 = x ^ x Q and -y 0 = y = y Q . When this 
two-dimensional exposure level passes through an imaging system which is relatively 
free from optical aberrations and distortions, maximum light intensity occurs at the 
center of the resulting image. That is, peak signal occurs where x Q and y Q approach 
zero. 


S rpte>t> s ) = fi( E b) ^Vj y ° f 0 P r 00 AE r (u,v;g,t) Gj f s (u,v;s) e i2 *( ux+v y) dv du dx dy ( 2 6) 

x o —0 ^ x 0 y 0 -y Q -x Q v -°o ’ ’ 


y 0 - 


-o 


Integrating with respect to x and y, and applying the limits, the peak signal becomes 


p OO p OO 

S r pfe,t, s ) = ^( Eb ) x liln 0 J ^ AE r (u,v;g,t) Gj >fjS (u,v;s) 


sin 27 tx 0 u sin 27ry 0 v 
27tx 0 u 27ry 0 v 


du dv 


y Q -o 


= AE ( Eb )I oo J oo AE r( u ’ v ;g- t ) G Z)f)S (u,v;s)du dv 


(27) 
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The peak signal value of the dark area may be found similarly and added to expres- 
sion (27) to yield a peak-to-peak signal 


Sp_ p (g,t,s) 



AE(u,v;g,t) G z>f>s (u,v;s)du dv 


( 28 ) 


where 

AE(u,v;g,t) = AE r (u,v;g,t) + AE q (u,v;g,t) 

The integral in the above equation presents a summation of the product of input 
signal frequency spectrum and system MTF. A film readout aperture scan across the 
film image transduces the two-dimensional distribution to a one- dimensional function. 

It is apparent that maximum signal will result if scanning occurs along sunline as illus- 
trated in figure 2(d). In pursuit of a simple to evaluate, approximate expression, a 
continuous square wave is considered rather than the isolated square-wave pulse. Of 
the resulting discrete frequency spectrum, the fundamental frequency component k is 
selected to estimate peak-to-peak signal values as indicated 

Sp_ p (g,t,s) = AE(g,t) G z>f>s (k;s)||(E b ) (29) 

Adding the noise term of equation (17), an approximate expression of the image signal- 
to-noise ratio becomes 

^£(g,t,s) = AE(g,t) ^^-^(Eb) (30) 

Change of exposure is related to the average bright and dark photometric value of 
the cone for a given exposure time (see eq. (24) and table I) by 

AE(g,t) = ^ A$(g) = 690tp A ¥(g) (31) 

4f n 2 

The quantity A$(g) is given in figure 5. The approximate cone model area (depicted 
in fig. 2(c)) is 

A i = 2 (f)(§) 2 = 0-39d 2 m 2 (32) 

The scanning aperture area referred to the lunar surface, using the scaling factor given 
previously, is 

Mi)(M) 2 = oW < 33 > 
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And the fundamental cone frequency is 

k = 1 line _ 76 lines /o 4 \ 

0.0132d mm d mm ' 

Finally, substituting equations (31) to (34) into equation (30), and considering that a peak- 
to-peak signal-to-noise ratio of 3 has been specified as detection threshold, a relation 
for the minimum detectable cone diameter can be written as 


* 


* 


\ 


* 


3 = [690tp A$(g)] 


Jo.39d 2 ^ 

rojL rj 





( 35 ) 


Grouping exposure- and phase-angle-dependent parameters on the left side of the equality 
sign and the frequency-dependent parameters on the right side, the equation becomes 


0.003 

tp A^(g) AT {— (E b ) 



(36) 


This equation has to be solved for d by trial and error. Solutions for exposure times 
of 1/25 and 1/50 second and a lunar albedo of 0.07 are plotted in figure 8. In one case, 
image smear is assumed to be absent; in the other case, smear rate is assumed to be a 
la value of 10 meters/second (ref. 16). 

Expe cted val ue of^ minimum detectable c one diameter .- The expected value of the 
minimum detectable cone diameter can be found approximately from the results given in 
figure 8 by using a numerical summation for equation (19) as follows: 

m n 

4*min) - II dmin (log Ei,gj) p(log Ei/gj) p(gjj (37) 

j=l i=l 

Exposure is replaced by log exposure for convenience as pointed out previously. 

The accuracy of these summations depends, of course, on the number of intervals 
used. (d m i n ^ is calculated here for n = 5 and m = 7 intervals ranging between the 
3a values of exposure and phase angle as depicted in figure 4. Results of these calcula- 
tions are shown in figure 9. 

Confidence levels in detection of a specifi ed cone diameter .- Confidence levels in 
the detection of specified cone diameters can be found by using equation (23) and the 
results of equation (37). Again, it is easier to work with log exposure instead of expo- 
sure. Results of these calculations for a specified 2-meter-diameter cone at no image 
smear and at a smear rate of 10 meters/second are shown in figures 10(a) and 10(b), 
respectively. Assuming that a shutter speed of 1/25 second is used for g > 74° and of 
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1/50 second for g < 74°, the confidence levels in detecting 2-, 3-, and 4-meter- 
diameter cones at the image smear rate of 10 meters/second are plotted in figure 11. 

Discussion of Results 

Several factors of primary importance may be noted from the results of these 
calculations. 

(1) As shown in figure 8(a), best picture quality in revealing small lunar detail 

occurs at high phase angles (70° to 83° for a lunar albedo of 0.07) and long exposure 
times (1/25 second) if image smear is negligible. Furthermore, picture quality degrades 
rapidly at phase angles above 83° and slowly below 70°. In the presence of some image 
smear, as shown in figure 8(b), shorter exposure times (1/50 second) and, therefore, * 

lower phase angles (62° to 74°) are preferable. 

(2) Curves depicting the variation of minimum detectable cone diameters with 
changes in phase angle (fig. 8) and those depicting its average value (fig. 9) differ in 
proportion to the rate of degradation with changes in exposure and phase angle and the 
magnitude of uncertainties in exposure and phase angle. 

(3) While a cone of 2-meter-diameter may be detected in a range of phase angles 
extending from 51° to 86° (fig. 8(a)) for a vehicle perilune of 46 kilometers, in the 
absence of image smear, and under expected conditions in lunar reflection and vehicle 
performance, confidence levels in the detection of cones this size decrease rapidly below 
0.5 outside the range of 54° to 85° (fig. 10(a)). 

(4) Final results are summarized in figure 11. Shown are the variations of con- 
fidence levels in the detection of 2-, 3-, and 4-meter-diameter cones with phase angle. 
Trade-offs in phase angle (trajectory) design are readily apparent between low confi- 
dence levels in the detection of small detail and higher and wider confidence bands in the 
detection of larger detail. 


CONCLUDING REMARKS 

An analysis is presented which relates photographic quality in revealing small 
lunar detail to viewing geometry and film exposure and hence to mission controllable 
parameters. Results of this analysis show that the best performance can be achieved 
for Lunar Orbiter’s vertical photography at high phase angles and slow shutter speeds 
providing that no error exists in the lunar photometric model and no deviation occurs 
from the expected vehicle operation. 

To consider uncertainties in lunar reflectance and system operation, a strategy is 
developed which permits the selection of vehicle trajectories and camera shutter speeds 
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for optimum confidence levels in achieving a specified performance. Application of this 
strategy leads to the selection of a range of phase angles extending from 65° to 80° as 
giving the highest level of confidence in the detection of 2-meter-diameter cones. A low 
and narrow confidence band in the detection of small lunar detail can be traded for a 
higher and wider confidence band in the detection of larger detail by small vehicle tra- 
jectory changes. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 26, 1967, 

814-11-00-03-23. 
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APPENDIX 


CALCULATION OF AN AVERAGE PHOTOMETRIC VALUE FOR THE 
BRIGHT AND DARK SIDE OF A LUNAR CONE 


The average photometric function of the bright and dark halves of a lunar cone 
depends on phase angle g and an effective lunar slope /3 e . This slope is measured 
around the right- circular cone as a function of the angle 0 away from the plane formed 
by the incident sun and the camera as depicted in figures 2 and 12 and may be written in 
terms of the coordinate variables as 


0 e (0) = -tan 


-II 




The equation for a right- circular cone is 


x 2 + y 


2 _ z 2 

n2 


(Al) 


(A2) 


where c is a constant. Solving this equation for z and differentiating with respect 
to y yields 

a/2 


z = c (x 2 + y 2 ) 


dz _ cy 
dy ~ i\V% 


(x 2 + y2) J 


Substituting x = sin 0 and y = cos 0 into equation (A4) gives 

dz 

dy 


^ = c cos 6 


When 0 = 0, 


dz _ 2h 
dy d 


where h is the cone height and d is the base diameter. Thus, 


dz 2h Q 
— = — — cos 6 
dy d 


The effective slope is, therefore. 


/3 = -tan" 1 ^ cos 0 
e |_d 


(A3) 

(A4) 


(A 5) 


(A6) 


(A7) 
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For a 26.6° cone, which has a base-to-height ratio of 4:1, this becomes 

j3 e = -tan”^[j}.5 cos 0] (A8) 

The average photometric function- of the bright half of the cone 4> r (g) and the dark 
half 4>q(g) may now be obtained from the photometric curves as indicated 


n 

<Mg) = £ ^ $ ife’£e) 
i=0 

(-90° < 9i < 90°) 

(A9a) 

n 

$q(g) = \ ^ $ife>/3 e ) 
i=0 

(90° < 0i < 270°) 

(A9b) 


Results of this averaging for n = 19 (10° intervals) is shown in figure 5. 
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TABLE I.- NOMINAL VALUES AND TOLERANCES 


Functions 

Nominal values 

3ct exposure 

3 a log exposure 

Solar constant, K 

1.3 x 10 5 meter-candles 


— 

Static lens transmission, r 

0.66 


— 

f- number, f n 

5.6 


— 

Phase angle, g 

0° - 90°, (±3.5°, 3cr) 


— 

Film 

See figure 9 

* ±23% 

0.10 

Shutter time (static), t 

1/25 and 1/50 sec 

±25% 

.11 

Shutter time (dynamic), t 

1/25 and 1/50 sec 

±5% 

.02 

Lunar albedo, p 

0.05 - 0.13 

±40% 

.17 

Photometric function, $ 

See figure 3 

±30% 

.13 

Exposure, E 

rss 

±64% 

0.28 
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Figure L- Block diagram of Lunar Orbiter photographic system. 
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(b) Phase angle probability distribution, 3 a = ±3.5°. 

Figure 4.- Gaussian log exposure and phase angle probability distribution approximated over discrete intervals 
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(a) Modulation transfer function for lens and film. 
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Modulation transfer function for image smear. 
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(c) Combined modulation transfer function for lens, film, and image smear. 


Figure 6.- Modulation transfer functions for photographic subsystem and image smear. Smear is assumed Gaussian and is designated by its 
standard deviation referred to the film image and to the lunar surface for a vehicle perilune of 46 km. 
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(a) Image smear absent. 



30 40 50 60 70 80 90 

Phase angle, deg 


(b) Image smear rate, 10 meters/second. 

Figure 9.- Variation of expected value of the minimum detectable 26.6° cone diameter with phase angle. Image smear is assumed Gaussian and is 
designated by its standard deviation referred to the lunar surface for a vehicle perilune of 46 km. Lunar albedo, 0.07. 
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Confidence levels Confidence levels 



(a) Image smear absent. 



Phase angle , deg 


(b) Image smear rate, 10 meters/second. 

Figure 10.- Variation of confidence levels in the detection of 2-meter-diameter, 26.6° cones with phase angle. Image smear is assumed 
Gaussian and is designated by its standard deviation referred to the lunar surface for a vehicle perilune of 46 km. Lunar albedo, 0.07. 
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